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In this paper, the Dvah-Gabadadze-Porrati (DGP) brane model is confronted by current cosmic 
observational data sets from geometrical and dynamical perspectives. On the geometrical side, the 
recent released Union2 557 of type la supernovae (SN la), the baryon acoustic oscillation (BAO) 
from Sloan Digital Sky Survey and the Two Degree Galaxy Redshift Survey (transverse and radial 
to line-of-sight data points), the cosmic microwave background (CMB) measurement given by the 
seven-year Wilkinson Microwave Anisotropy Probe observations (shift parameters R, la{Zi,) and 
redshift at the last scatter surface z*), ages of high redshifts galaxies, i.e. the lookback time (LT) 
and the high redshift Gamma Ray Bursts (GRBs) are used. On the dynamical side, data points 
about the growth function (GF) of matter linear perturbations are used. Using the same data 
sets combination, we also constrain the flat ACDM model as a comparison. The results show 
that current geometrical and dynamical observational data sets much favor flat ACDM model and 
the departure from it is above 4(t(6(t) for spatially flat DGP model with(without) SN systematic 
errors. The consistence of growth function data points is checked in terms of relative departure of 
redshift-distance relation. 

PACS numbers: 98.80.Es, 04.50.-h, 04.50. Kd, 95.36.+X 



I. INTRODUCTION 

Understanding the current accelerated expansion of our universe has become one of the most important issues of 
modern cosmology P, However, as so far, we still know little about the nature of current accelerated expansion. 
In general, from the phenomenological points of view, the possible models can be classified by the form of Friedmann 
equation. One is 

- /(p), (1) 

where the extra energy component(s) is(are) added on the right hand of Einstein's equations, and f{p) is a function 
of energy density p which can be composed from dark matter and extra energy components. Here H is the Hubble 
parameter. The other is 

9iH') = p, (2) 

where the energy component (s) is (are) invariant and the Gravity Theory or equivalently the relation of Hubble 
parameter H with the conventional matter (dark matter at late time) can be altered, g{H^) is a function of or 
H. For recent reviews about dark energy and an accelerated expansion universe, please see [3|. 

In this paper, we will consider a leading modified gravity model, Dvali-Gabadadze-Porrati (DGP) brane model 
for the review, please see [5]. In DGP brane model, a tensionless four-dimensional brane (a hypersurface with a 
vanishing cosmological constant) is embedded in a five dimensional bulk which is a flat Minkowski space-time, where 
the gauge forces are confined on the brane and gravity can propagate in all dimensions freely. Below the crossover scale 
Tc, the gravity appears four-dimensional. However, above the scale r^ the gravity can leak into the extra dimension 
and make the conventional four-dimensional gravity altered. It is due to the leakage of gravity, the current universe 
appears an accelerated expansion phase. 

In the past years, the DGP model has been constrained by cosmic observational data sets, for recent constrained 
results, please see [a] in geometrical side and in dynamical side. Importantly, Lombriser et.al. Q constrained DGP 
brane model exhaustively by adopting a parameterized post- Friedmann description of gravity, where all of the CMB 
data, including the largest scales, and its correlation with galaxies in addition to the geometrical constraints from 
supernovae distances and the Hubble constant were utilized. Of course, it is more important and necessary to constrain 
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a cosmological model by using all of CMB data and its correlation with large scale structure, say galaxy. However, it is 
much complicated when we have to resort to some kinds of Boltzmann-codes, say the famous CMBfast and CAMB. 
It would be interesting, if this complication is avoided, and a relative stronger prediction is arrived. In this paper, we 
will investigate this kind of possibility by combining the geometrical and dynamical perspectives without the above 
mentioned complication. Via this combination, the constraint becomes much tight and efficient because it relies on 
both gravity theory and background geometry. On the geometrical side, we will consider the luminosity-distance 
relation of SN and GRB, the standard rulers from BAO, x-ray gas fraction and CMB, also the lookback time of high 
redshift galaxy. For the dynamical perspective, we mainly consider the growth function S{z) = ^(z) of the linear 
matter density contrast as a function of redshift. As well known, the data points of growth function were obtained 
with the help of ACDM model and the reliability will be reduced when it is used to constrain other cosmological 
models, however this problem can be evaded in the DGP model if it has similar expansion history to that of ACDM 
0, Hj. In fact, the use of redshift-distance relation of ACDM (with flm — 0.25 for the data point at z = 0.77 and 
rim = 0.30 for the remained five data points) makes the the data points of growth function weak to constrain other 
cosmological model In @, to use these data points to find the growth index of DGP model, the DGP model 

parameters were selected properly to make redshift-distance relation be the same as that of ACDM model. In fact, we 
can firstly neglect this weak point and just use these data points as dynamical constraint beyond ACDM model. After 
it was done, we can check the consistence via investigating the possible departure from ACDM model with the best 
fit values of model parameters. If the departure is large, the data points must be discarded and the data fitting is not 
reliable. In fact, one will find that the redshift-distance relations of ACDM and DGP model are almost the same and 
the relative departure from ACDM model is up to 10% in terms of Hor(z) when ilm varies in the range [0.28,0.32], 
where flm = 0.30 is fixed in ACDM model. When the values of ilm increase, the relative departure becomes small. 
And the departure is smaller than the errors of growth function data points. Based on this point, the data points of 
growth function can be used safely to constrain DGP model. Also, when the data points of growth function are used 
to confront other cosmological models, the consistence must be checked. 
For the DGP model, the modified Friedmann equation is given as [llj 



3Af2 4r2 2rc 



(3) 



where Vc = M'^^/^M^ is the so-called crossover scale, Mpi and are four- and five-dimensional reduced Planck mass 
respectively, p is the cosmological fluid which includes conventional matter contents and radiation, fc = 0, ±1 is three 
dimensional spatial curvature factor. In terms of the dimensionless density parameters, the Friedmann equation can 
be rewritten as 

E'^{z) ^nk{i + zf + (^^^^(i + z)3 -t- nr{i + zf + rir^ + y/a^y , (4) 

where E'^(z) = H'^{z)/Hq and ilr^ = l/(4r^iJg) is a constant which respects to the constraint equation 



f7fe + [^^nrn + Or + ^r^ + =1- (5) 

For the spatially flat case (fife = 0), the above equation reduces to ilr^ = (1 ~ ~ r2r)^/4. 



II. METHOD AND RESULTS 

In our calculations, we have taken the total likelihood function L oc e"-*^^/^ to be the products of the separate 
likelihoods of SN (with and without systematic errors), BAO, CMB, GRBs, CBF, LT and GF. Then we get 

= xIn + X%AO + xImb + X%RBs + XcBF + Xlt + Xgf, (6) 

where the separate likelihoods of SN, BAO, CMB, GRBs, CBF, LT, GF and the current observational datasets used 
in this paper are shown in the Appendix [A] 

In our analysis, we perform a global fitting to determine the cosmological parameters using the Markov Chain 
Monte Carlo (MCMC) method. The MCMC method is based on the pubhcly avaflable CosmoMC package ^ and 
the modified CosmoMC package [T3l - [l5l |. including the X-ray cluster gas mass fraction. For our models we have 
modified these packages to add some subroutines, for example the code calculating the growth function likelihood 
etc. The following basic cosmological parameters (ftbh^, i^ch^) are varying with top-hat priors: the physical baryon 
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Datasets 


Parameters Xmin/d. o.f 




Ho 


GF5 

GF6 
GF5(SN Sys.) 
GF6(SN Sys.) 


9 697.342/645 
9 707.750/646 
9 663.561/645 
9 673.641/646 


U.^UO_f, Q3Q2 

U.ZOO_o Q3Q4 

U.ZOO_o 0286 
n 9Q7+0.0367 
U.Z»(_o 0391 


65.0873t?;4i3 
63.754l|^*^ 

U0.1Z1_2 0471 



TABLE I: The results of Xmini and Ho with la regions are Usted. GF5 denotes 5 growth function data points are used 
without z = 3.0 data point. And, GF6 denotes the corresponding 6 data points case. SN Sys. denotes the resuhs with SN 
systematic errors, d.o.f denotes the degrees of freedom. 



Datasets 


Parameters 


xLn/d.O.f 


^T7l 


Ho 




non-flat 
non-flat(SN Sys.) 


10 
10 


684.251/644 
647.543/644 


n O7n+0'0278 
U.Z(U_o.0323 
n 9Qc;+0.0419 
U.Z»0_o.o345 


66.798t^:f83 
65.519t|.«« 


n ni9q+o oo''89 

U.UlZO_o 00993 
U.U100_o. 01052 



TABLE IL The results of Xmim ^m, Hq and Qk without z = 3.0 data point in non-flat DGP model. SN Sys. denotes the 
results with SN systematic errors, d.o.f denotes the degrees of freedom. 



density ftbli^ £ [0.005,0.1], the physical dark matter energy density ftch^ e [0.01,0.99]. And, in the data fitting 
process another seven parameters (if, 77, 7, 6o, Q^b, so: Q^s) included in the X-ray gas mass fraction fgas are treated as 
free parameters. As a byproduct the best fitting values of these parameters are obtained. And, these values can also 
be taken accounted as a check of data fitting. 

From the observational results listed in Tab. IVIIl one has noticed that the datum at 2; = 3.0 is odd in some degree 
for its value is above 1. So in data fitting process, we have considered two cases with and without inclusion of z = 3.0 
data point. The constrained results are shown in Tab. Uand Fig. [T] (with SN systematic errors). From Tab. HI one 
can easily see that the value of Xmin smaller and the constraint is tighter without z = 3.0 data point than that 
with it. Under this status, one would constrain a cosmological model without the inclusion oi z = 3.0 data point. 
For the non-flat DGP model, we summarize the model parameters values in Tab. Ulland Fig. O (with SN systematic 
errors). As a comparison, using the same combination (with SN systematic errors) to the flat ACDM model, we have 
the results X™„ = 626.057, f7„ = 0.28ll^;!^i^^^ and Hq = 69.658l|?4^^ One can see that the best fit values of rj„ 
and la errors of model parameters are enlarged when SN systematic errors are included. As contrasts, the current 
Hubble values and Xmin lower in the case of SN systematic errors included. 



FIG. 1: The 2-D contours with la, 2a regions and 1-D marginalized distribution of Q,m, Ho- Left Panel: Xmin ~ 663.561, 
Vim ~ 0.285lo'o286 ^'id Ho ~ 63.75412 575 with 5 growth function data points, z = 3.0 is removed. Right Panel: Xmin ~ 673.641, 
= 0.297io;o39l and Ho = 63.12112:^4^^ with aU 6 growth function data points included. 

Now, it's time to check the consistence of our treatment of the points of growth function. We calculate the relative 
departure of comoving distance 



rel.dep. = 

rhCDMKz) 



(7) 
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FIG. 2: The 2-D contours with la, 2a regions and 1-D marginalized distribution of Qrn, Ho and flk in non-flat DGP model 
without z = 3.0 data point. 

where ftm — 0.30 is fixed in ACDM model. The result is shown in Fig. |3l In fact, one can check the consistence via 
the relative departure in terms HQr{z)/c with the variables of z and fim, the corresponding result is shown in a 3D 
plot, see Fig. El Clearly, the relative departures are up to 10% for neglecting the uncertainty of Hq. And, with the 
increasing value of ftm, the departure is amplifying. In fact, as shown in Fig. |31 the departure is just up to 8% when 
the best fit values of model parameters are adopted. 
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0.00 I ' ' ' ' ' ' ' ' ' < ' ' ' ' ' ' ' ' ' 1 
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FIG. 3: The relative departure of redshift-distance r{z) of DGP model with the best fit values of model parameters from ACDM 
model with fim = 0.30 and Ho = 69.658. 
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III. CONCLUSION AND DISCUSSION 



In summary, in this paper we have performed a gfobal fitting on the cosmological parameters in both the flat DGP 
model and the non-flat DGP model by using a completely consistent analysis from the geometrical and dynamical 
perspectives. On the geometrical side, the X-ray gas mass fraction observation, type la supernovae data from Union2 
set, transverse and radial baryon acoustic oscillations data from SDSS, the measurement data on current Cosmic 
Microwave Background from the seven-year WMAP observations, the lookback time data derived from the ages of 
galaxy and clusters and the high redshift gamma ray bursts. On the dynamical side, the data points about the 
growth function of matter linear perturbations are included. The constrained results are shown in Tab. |T] for the 
flat case and Tab. |TT] for the non-flat case. The results show that the z — 3.0 data point of growth function 
is odd in some degrees. So, in the future work when the observational growth function data points are used as 
cosmic constraint, this point would be removed. As a comparison, using the same data points combination, the 
flat ACDM model was constrained, where we have Xm„ = 626.057, Q.^ = 0.28ll^;[5i^^^ and = 69.658t^;?4^^. 
Clearly, one can easily find that A^^ — 37.504 which in terms of a distance is 4.25 for 9 parameters via the formula 
1 — r(z^/2, /S.-}^ /2) /T{i> /2) = Erf ((icr/"\/2), where = 9 is the number of free model parameters. It means that current 
geometrical and dynamical observational data sets favor flat ACDM model more above 4cr than that for spatially flat 
DGP model. This conclusion is consistent and comparable with that of |8|. When the same process of data fitting is 
implemented without SN systematic errors, the departure from ACDM model can be improved to about 6tT for the 
flat DGP model. To check the consistence of growth function data points, the relative departure of DGP model with 
best fit model parameters from ACDM model with fixed Vl — 0.30 is up to 8%. In general, the relative departure in 
terms of HQr{z) is up to 10% when varies in the range [0.28, 0.32] and fim — 0.30 is fixed in ACDM model. With 
this observation, we can say that growth function data points can be used to constrain DGP model. So far, we have 
seen that it is possible to give some strong prediction without having to resort to the complicated modification of the 
Boltzmann-codes. At least, for flat DGP model, it is possible. It is because that the discrepancy will be enlarged 
when more observational data points are included. We expect this work can shed light on distinguish dark energy 
models. At last, we are pleasant to give some comments on the growth function data points used to constrain other 
cosmological models. The key point is the dependence of distance-redshift relation of spatially flat ACDM model with 
fixed value of flm — 0.30. Of course, one can firstly neglect this fact and do data fitting as what we have done in 
this paper. Then, one has to check the consistence. In fact, we can get around this problem in a different way. It is 
the full consideration of distance-redshift relation and introduction of some kinds of uncertainties in terms of redshift 
and growth function data points, put in another words, via the possible introduction of extra errors on redshifts and 
growth function data points. Of course, it is out the range of this paper. In the future work, this possibility will be 
researched. 
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Appendix A: Cosmological Constraints Methods 

1. Type la Supernovae constraints 

Recently, SCP (Supernova Cosmology Project) collaboration released their Union2 dataset which consists of 557 
SN la flG*!. The distance modulus fi{z) is defined as 

fith{z)^5\ogj^„[dL{z)]+ ^iQ, (Al) 

where ^^(z) is the Hubble-free luminosity distance Hi^dL{z) / c — HQdA{z){l + z)'^/c, with Hq the Hubble constant, 
defined through the re-normalized quantity h as Hq — 100ft- km s^^Mpc^ , and /iq = 42.38 — 51og]^g h. Where d^lz) 
is defined as 



di(z) = (1 + z)r(z), r(z) = — ^==sinn 



dz' 



(A2) 



/o E{z')_ 

where E'^{z) = H'^{z)/Hq. Additionally, the observed distance moduli jiobsizi) of SN la at Zi is 

tJ^obs[zi) = mobs[zi) - M, (A3) 

where M is their absolute magnitudes. 

For the SN la dataset, the best fit values of the parameters Ps can be determined by a likelihood analysis, based 
on the calculation of 

X^{p M') = ^ {l^obs{zi) ~ lJ.th{Ps,z,)f 
SN 

E {51ogio[rfL(Ps,^:»)] -mobsjzi) +M'y 
SN ^ 

where M' = /^o + M is a nuisance parameter which includes the absolute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analytically [l7j as 



/ + 00 
exp 
CO 



-f-oo 

' dM\ 



resulting to 



with 



{51ogio[dL(Ps, Zj)] - moi,s{zi)y 

SN — 

^ _ ^ 5\ogi„[d L{ps, Zi)] - mobs{zi) 



SN 
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Relation (|A4[) has a minimum at the nuisance parameter value M' = B/C, which contains information of the values 
of h and M . Therefore, one can extract the values of h and M provided the knowledge of one of them. Finally, note 
that the expression 

xIj,{Ps,B/C)^A-{B^IC)^ (A7) 

which coincides to Eq. (jASp up to a constant, is often used in the likelihood analysis [13; [13, and thus in this case 
the results will not be affected by a flat M' distribution. It worths noting that the results will be altered without 
the systematic errors. In this work, two cases with and without systematic errors are considered together. When the 
systematic errors are included, the corresponding A, B, C are expressed as 

SN 

^ = ^ {51ogio[<iL(Ps,Zi)] - mobs{zi)} ■ Cj^/ ■ {5 log j^Q[dL{ps,Zj)] - niobsizj)} , 

i,j 
SN 

i 

SN 

C = (A8) 

i 

where is the inverse of covariance matrix with systematic errors. For the details and covariance matrice, one can 
find them in Ref. IJJ^ and the web site ^, where one can also find the covariance matrix without systematic errors. 
Our form (jA6l) is equivalent to (jA8l) when is the inverse of covariance matrix without systematic errors. 



2. Baryon Acoustic Oscillation constraints 

The Baryon Acoustic Oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, and measure the distance-redshift relation at z = 0.2. Additionally, Baryon Acoustic Oscillations in the 
clustering of the SDSS luminous red galaxies measure the distance-redshift relation at z = 0.35. The observed scale 
of the BAO calculated from these samples, as well as from the combined sample, are jointly analyzed using estimates 
of the correlated errors to constrain the form of the distance measure Dv{z) |19l - l22j 

dz' 



The peak positions of the BAO depend on the ratio of Dv{z) to the sound horizon size at the drag epoch (where 
baryons were released from photons) Zd, which can be obtained by using a fitting formula ,23j|: 

with 

hi = 0.313(a„7i')"°-^'^[l + 0.607(an/i')°-^'^], (All) 
&2 = 0.238(17™/i^)°-223. (A12) 

In this paper, we use the data of rs{zd)/Dv{z), which are listed in Table Hill where rs{z) is the comoving sound 
horizon size 



/■* Csdt f Csda , 
rs[z) = c I — = c I -Y-rp = c / dz 



c 



Q a Jq a?H H{z) 
1/(1+^) da 



VsJq a^H{a)./lTWhJWh)a}' 



(A13) 



^ http:/ /supernova. lbl.gov/Union/ 
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where Cg is the sound speed of the photon— baryon fluid [2J-|26j : 



and here = 2.469 x lO^^/i"^ for Tcmb = 2.725if 



z 


rs{zd)/Dv{z) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE III: The observational rs(zd)/Dv{z) data [23]. 
Using the data of BAO in Table Hill and the inverse covariance matrix V^-^ in (20j : 

T/-=f 30124.1 -17226.9 y 

\^ -17226.9 86976.6 J ^ ^ 

The radial (line-of-sight) BAO scale measurement from galaxy power spectra give constraint to cosmological pa- 
rameters via the relation 

AzsAoiz) = ^^^^ (A16) 
c 

at two redshifts z = 0.24 and z = 0.43, the corresponding values are Azbao{z — 0.24) — 0.0407 ± 0.0011 and 
AzBAo{z = 0.43) = 0.0442 ± 0.0015 respectively 
Thus, the X%ao(Ps) is given as 

xlAoiPs) - X^V-^X + \^^BAoiz^0.24) -0M07r ^ [A...o(z = 0.43) - 0.0442]^ 

A^BAOyP^I -r- 0.00112 ^ 0.0015^ ' ^ ' 

where AT is a column vector formed from the values of theory minus the corresponding observational data, with 

/ _ 0.1905 ^ 



\ Dv(0.35 

and A* denotes its transpose. 



X = ^//o^?) , (A18) 

\ r,{za) - 0.1097 / ^ ^ 



3. Cosmic Microwave Background constraints 

The CMB shift parameter R is provided by [H] 

i?(z,) = V^^sinn[/^ -^], (A19) 

here, the redshift (the decoupling epoch of photons) is obtained by using the fitting function [2^ 

2* = 1048 [1 + 0.00124(r2fa/i2)-" '^3S] [l + ^ 
where the functions gi and ^2 read 

51 = 0.0783(17b/l2)-0.238 (^^ 39 5(^^^^2)0.763)-!^ 

52 = 0.560 (l± 21. 1(17b/i2)i «i)"\ 

In addition, the acoustic scale is related to the distance ratio and is expressed as 



dz' 
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7 - year ML 


7 — year mean 


error, a 




302.09 


302.69 


0.76 




1.725 


1.726 


0.018 


2, 


1091.3 


1091.36 


0.91 



Using the data of Ia,R,z^ in [sS], which are hsted in Table ITVl and their covariance matrix of [Ia{z^,), R{z^), z^] 
referring to [30i] : 



2.305 29.698 



1.333 



C-^ = I 29.698 6825.270 -113.180 | , 
-1.333 -113.180 3.414 



we can calculate the likelihood L as Xcmb ~ ~ 

21nL: 



(A21) 



(A22) 



where Ad^ = di — df is a row vector, and di — {I a, R, z^). 



4. Gamma Ray Bursts 

Following [sTl, we consider the well-known Amati's Ep^i—Eiso correlation [s^ - lssj in GRBs, where Ep^i = Ep^obs{'^+z) 
is the cosmological rest-frame spectral peak energy, and Eiso is the isotropic energy 



Ei, 



ATrdlS,oio/il + z) 



(A23) 



in which dj^ and S\,oio £^re the luminosity distance and the bolometric fluence of the GRBs respectively. Following 
l3ll. we rewrite the Amati's relation as 



1 Eiqq Ep -i 

log = a + b log ■ 



erg ' 300keV 

In [36j . Wang defined a set of model-independent distance measurements {fp(zi)}: 

rp(zo) z c 



(A24) 



(A25) 



where r(z) = di^{z)/{l + z) is the comoving distance at redshift z, z^ is the lowest GRBs redshift. Then, the 
cosmological model can be constrained by GRBs via 



Xgr.Bs{Ps) = [Arp(z,)] • (Cowg^j^),^ • [Arp(zj)], 
^rp{zi) = ff*'^(z,)-fp(z,). 



(A26) 
(A27) 



where fp{zi) is defined by Eq. (IA25P and {Cov^^g)ij , i,j = l...A^ is the covariance matrix. In this way, the constraints 
from amount observational GRBs data are projected into the relative few quantities fp{zi),i = 1...N. 

Following the method proposed by Wang [36J, Xu obtained = 5 model-independent distances data points and 
their covariance matrix by using 109 GRBs via Amati's Ep^i — Eiso correlation [37|. The resulted model- independent 
distances and covariance matrix from 109 GRBs are shown below in Tab. IVland Eq. (|A29|) . The {fp{zi)}, i — 1, 5 
correlation matrix is given by 



[Covgrb) 



( 1.0000 0.7780 0.8095 0.6777 0.4661 \ 

0.7780 1.0000 0.7260 0.6712 0.3880 

0.8095 0.7260 1.0000 0.6046 0.5032 

0.6777 0.6712 0.6046 1.0000 0.1557 

V 0.4661 0.3880 0.5032 0.1557 1.0000 / 



(A28) 
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z 




a{f,{z)Y 


CT(r-p(2)) 





0.0331 


1.0000 


— 


— 


1 


1.0000 


0.9320 


0.1711 


0.1720 


2 


2.0700 


0.9180 


0.1720 


0.1718 


3 


3.0000 


0.7795 


0.1630 


0.1629 


4 


4.0480 


0.7652 


0.1936 


0.1939 


5 


8.1000 


1.1475 


0.4297 


0.4389 



TABLE V: Distances measured form 109 GRBs via Amati's correlation with Icr upper and lower uncertainties [37 
as lowest redshift was adopted. 



20 = 0.0331 



and the covariance matrix is given by 



[CovGRB)ii = o{rp[zC))a[ry{zj)){CovGRB)i3, 



where 



a(r^{zO) ^ a{f^{z,))- , if fp(z) < fp(z)'^"*^ 
the cr(fp(zi))+ and (7(fp(zi))~ are the Itr errors Usted in Tab. IVl 



(A29) 



(A30) 
(A31) 



5. The X-ray gas mass fraction constraints 



According to the X-ray cluster gas mass fraction observation, the baryon mass fraction in clusters of galaxies (CBF) 
can be utilized to constrain cosmological parameters. The X-ray gas mass fraction, /gas, is defined as the ratio of 
the X-ray gas mass to the total mass of a cluster, which is approximately independent on the redshift for the hot 
(fcT > 5fceT^), dynamically relaxed clusters at the radii larger than the innermost core ^2500- The X-ray gas mass 
fraction, /gas, can be derived from the observed X-ray surface brightness profile and the deprojected temperature 
profile of X-ray gas under the assumptions of spherical symmetry and hydrostatic equilibrium. Basing on these 
assumptions above, Allen et al. [l^ selected 42 hot (fcT > 5fceF), X-ray luminous, dynamically relaxed clusters for 
jgas measurements. The stringent restriction to the selected sample can not only reduce maximally the effect of the 
systematic scatter in /gas data, but also ensure that the /gas data is independent on temperature. In the framework 
of the ACDM reference cosmology, the X-ray gas mass fraction is presented as [l^ 



KA-ih{z) 
1 + s{z) 









dA{z) 



1.5 



(A32) 



where A is the angular correction factor, which is caused by the change in angle for the current test model ^2500 in 
comparison with that of the reference cosmology 6*2560*^^' 



A- 



/flACDM \ ^ 
I '^2500 > 



72500 



H{z)dA{z) 



[H{z)dA{z)]^^^^ 



(A33) 



here, the index 77 is the slope of the fgasir/r25oo) data within the radius r25oo, with the best-fit average value 
77 — 0.214 ± 0.022 13']. And the angular diameter distance is given by 



dA{z) = 



(1 + z)/^ 



:Sinn 



dz' 

H{^y 



(A34) 



where sinnn(A/|f7fe|a;) respectively denotes sin(-\/|r2fc|a;), y^TIfe|a;, sinh(-\/|f^fe|a;) for fi^ < 0, flfc = and fifc > 0. 

In equation (IA32P , the parameter 7 denotes permissible departures from the assumption of hydrostatic equilibrium, 
due to non-thermal pressure support; the bias factor b{z) — 6o(l + oib^) accounts for uncertainties in the cluster 
depletion factor; s{z) = So(l -I- asz) accounts for uncertainties of the baryonic mass fraction in stars and a Gaussian 
prior for sq is employed, with sq = (0.16 ± 0.()b)h^Q [l^; the factor K is used to describe the combined effects of the 
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residual uncertainties, such as the instrumental calibration and certain X-ray modelling issues, and a Gaussian prior 
for the 'calibration' factor is considered by K — 1.0 ± 0.1 ^131] : 

Following the method in Ref. [l3, H^l and adopting the updated 42 observational fgas data in Ref. 13], the best 
fit values of the model parameters for the X-ray gas mass fraction analysis are determined by minimizing. 



2 _ [fgas i^i) f gasi^^i)] 



x'cBF = > - , (A35) 



where ct/ ^S^i) the statistical uncertainties (Table 3 of [131). As pointed out in [l3|, the acquiescent systematic 
uncertainties have been considered according to the parameters i.e. 77, s{z) and K . 



6. Lookback Time 

Since the seminal work of Sandage [s^ who defines the lookback time as the difference between the present age of 
the Universe (to) and its age {tz) when a particular light ray at redshift z was emitted, one has used the lookback 
time-redshift relation to constrain cosmological models, for example lookback time as a constraint to DGP model (ioj . 
The lookback time-redshift relation is given by 



Following Capozziello et al. [4l[, one can define the age t{zi) of an object (e.g., a galaxy, a quasar or a galaxy cluster) 

ween 



at redshift Zi as the difference between the age of the Universe at Zi and the age zp when the object was born 

dz' r dz 



{l + z')H{z') J,^ {l + z')H{z') 
- tL{zF)-tL{z^)■ (A37) 
Then the observed lookback time to an object at Zi can be defined as 

'(^.) = tdzp) ~ t{z{) = [t°^' ~ t{zi)] - - tLizp)] 

= tf' -t{z,)-df, (A38) 

where df is the delay factor which accounts for our ignorance about the absolute age of universe when the object is 
formed at t{zF)- To constrain a cosmological model by using lookback time of an object, one can use 

2 _sr^ [tL{z,)-tt{z,,df)f [to-tg"-]^ 

Xage — 2 + ^2 (.AOaj 

•^T '-'tabs 

where a'^ = af + <y%bs , is the uncertainty of the individual lookback time to the i*'' galaxy of our sample and a^ots 
is the uncertainty on the total expansion age of the universe. After marginalizing the 'nuisance' parameter c?/, one 
can use the following method to constrain a cosmological model by using lookback time [42j 



/•oo 

xIt{Ps) = -2 In/ d(d/)exp(-x^,e/2) 



= A-^ + D - 21n[y^erfc(-^)], (A40) 



where 

^-E^^ ^ = E^^ (A41) 

where A is 

A = ti(z,)-[t°''^-t(z,)] (A42) 

and D is the second term of Eq. (jA39|, erfc(a;) = 1 — erf(x) is the complementary error function of the variable x. 
The observational data points of the age of galaxies are shown in Tab. IVII . The current observational universe age 
is t"^" = 13.75 ± 0.13Gyr [s^. 
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Zi 


U{zi) (Gyr) 


0.1171 


10.2 


0.1174 


10.0 


0.2220 


9.0 


0.2311 


9.0 


0.3559 


7.6 


0.4520 


6.8 


0.5750 


7.0 


0.6440 


6.0 


0.6760 


6.0 


0.8330 


6.0 


0.8360 


5.8 


0.9220 


5.5 


1.179 


4.6 


1.222 


3.5 


1.224 


4.3 


1.225 


3.5 


1.226 


3.5 


1.340 


3.4 


1.380 


3.5 


1.383 


3.5 


1.396 


3.6 


1.430 


3.2 


1.450 


3.2 


1.488 


3.0 


1.490 


3.6 


1.493 


3.2 


1.510 


2.8 


1.550 


3.0 


1.576 


2.5 


1.642 


3.0 


1.725 


2.6 


1.845 


2.5 


0.60 


9.20 


0.70 


9.80 


0.80 


3.41 


0.10 


3.08 


0.25 


4.84 


1.27 


12.13 



TABLE VI: Galaxy ages (author?) 'i^. The last 6 data points are taken from the Capozziello et al. [ZJ] (Table 1) ages of 6 
galaxy clusters in the redshift range 0.10 < z < 1.27. The Ict error are taken as 12% of the age of galaxy with exception of the 
last 6 data points. The Icr error of the last 6 data points are set as 1. 



7. Growth Function of Matter Linear Perturbations 



The linear growth factor of the matter density perturbation D{cl) is defined as 



D{a) = , " , (A43) 
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which is subject to the evolution equation [38|, |4J| 

'3 E'{a) 



D {k,a) 



E{a) 



D'ik,a) - ^-^^f(k,a)Dik,a) = 



(A44) 



with initial conditions D{a) ~ a for a ~ on sub-Hubble scales, where ' denotes derivative with respect to scale factor 
a and E{a) — H{a)/HQ. The forms of f{k,a) depend on the dynamical equations of particular gravity theory. For 
general relativity, /(fc, a) = 1. For the case of DGP model, f{k, a) takes the form 



/(fc,a) = (l + — ), 



where a is 



a = 1 



a E'(a) 

1 H — 

3 E{a) 



The observational growth rate relates to D{a) via 



_ aD'{a) _ d\nS 



D{a) dlna 

The current data sets about growth function are listed in Tab. IVIIl 



o6s 



9i 



Reference 



0.15 
0.35 
0.55 
0.77 
1.4 
3.0 



0.51 ±0.11 
0.70 ±0.18 
0.75 ±0.18 
0.91 ±0.36 
0.90 ±0.24 
1.46 ±0.29 



[45] 
[46] 
[42 
[48] 

m 

[50] 



(A45) 



(A46) 



(A47) 



TABLE VII: Observed perturbation growth as a function of redshift z, where g{a) — 
Then, the constraint from growth function is in the form 



in terms of a. 



(A48) 
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